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^ ! Abstract 






We formulate a non-relativistic Hamiltonian in order to describe the interaction between a moving 
dielectric membrane and radiation pressure. Such a Hamiltonian is derived without making use 
of the single-mode adiabatic approximation and linear approximation, and hence it enables us to 

Qh| incorporate multi-mode effects and general (non-relativistic) motion of the membrane in cavity 

-I— ' ■ 

Ch ' optomechanics. By performing second quantization, we show how a set of generalized Fock states 






can be constructed to represent quantum states of the membrane and cavity field. In addition, we 
discuss examples showing how photon scattering among different cavity modes would modify the 



(N 

K»" ' interaction strengths and the mechanical frequency of the membrane. 
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I. INTRODUCTION 



The interplay between optical and mechanical degrees of freedom via radiation pressure 



is at the heart of cavity optomechanics 






M 



^ . With the recent advances in cooling techniques 



in optomechanical setups |5l-ll0l|. it is becoming possible to access quantum ground states, 



and study the interplay at the quantum level experimentally. This may 



plications in quantum information based on optomechanical coupling llNl5|. Specifically, 



ead to novel ap- 



the system formed by a movable dielectric membrane inside a optical cavity provides a basic 



configuration to explore quantum phenomena in macroscopic o 
classical mechanical states 



23. 



23| and cavity QED effects 
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mm 



jjects |l6l-l2lj. such as non- 



251. Technically, the system 



has the advantage that it enables a strong and tunable optomechanical coupling, with the 
possibility of achieving nonlinearity in the membrane's displacement [2l|. 

The theoretical analysis of a moving membrane system is mainly based on the Hamilto- 
nian (l9|-|25| 



^ I ^ dujix) 



dx 



a'a 



+ nnrnb^b, 



(1) 



^0, 



where a and b are the annihilation operators for optical and mechanical modes, respectively, 
Xm oc {b + b^) is the displacement operator of dielectric motion, and Qm is the natural 
frequency of the membrane. The radiation pressure coupling between the two degrees of 
freedom is contained in the membrane-position-dependent field frequency uj{xm), which can 
be linearized around an equilibrium position xq. If the first derivative of field frequency 
vanishes, then the coupling is dominantly quadratic in membrane displacements 
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m- 



We note that the essential assumption of the model (1) is that all photons would stay in 
the same cavity mode throughout the evolution because of the slow motion of the dielec- 
tric membrane. This is an adiabatic single-mode approximation in which photon scattering 
among different cavity modes are neglected. However, we point out that the mode cou- 
pling induced by membrane's motion can become significant in non-adiabatic regimes. This 
happens especially when the oscillation frequency of the dielectric is close to the frequency 
spacing between two cavity modes, in which case transitions between the two modes can 
be resonantly enhanced. Indeed, by exploiting such a kind of resonance, Dobrindt and Kip- 
penberg have recently indicated an optomechanical displacement transducer with a high 



sensitivity [26| . Hence a natural question of the moving-membrane system is how a Hamil- 
tonian model can be rigorously formulated, without employing the adiabatic single-mode 
approximation. Such a Hamiltonian would provide us with a basis of studying the quantum 
mechanics of field-membrane systems. A better understanding of the field-membrane inter- 
action also opens up possibilities of new schemes to manipulate quantum states of both the 
light field and mechanical motion. 

A key to address the field-membrane interaction is to treat both the field and the moving 
membrane consistently as dynamical variables. A similar problem for a moving perfect- 
mirror system has been treated in Ref . [27| . The case of moving-dielectric system has yet 
been formulated, although the field Hamiltonian for a dielectric or partially transparent 
mirror moving in a prescribed trajectory has been discussed in the context of dynamical 



Casimir effect 28|-|30|. A major conceptual difficulty of the problem is that the normal 
modes associated with cavity field depend on the position of the dielectric, which enters 
as a dynamical variable. These field modes change with time as the dielectric membrane 
moves, which in turn affect the radiation pressure on the membrane. A consistent approach 
to the coupled field-membrane dynamics is hence essential to tackle the problem. Recently, 
Biancofiore et al. [3l| have constructed a Hamiltonian based on a linearized form of the 
radiation pressure coupling (i.e., to first order in Xm) ■ While their Hamiltonian may also 
address non-adiabatic photon scattering among cavity modes, it remains unclear how the 
constructed Hamiltonian can be extended beyond the linear coupling. 

The main purpose of this paper is to provide a Hamiltonian formulation that describes 
the optomechanical coupling without using the single-mode adiabatic approximation and 
linear approximation. Based on the interaction between macroscopic dielectric and electro- 
magnetic field, we first derive the Lagrangian and Hamiltonian of the classical counterpart 
of the system. The canonical quantization of the Hamiltonian and the membrane-position 
dependent Fock states are then introduced. Our Hamiltonian indicates how photons de- 
fined by such Fock states can be coupled to various cavity modes through the motion of the 
mirror. In the regime where adiabatic single-mode approximation and linear approximation 
are applicable, our Hamiltonian can be reduced to the usual form (1). Near the end of this 
paper we indicate some physical consequences arising from the involvement of the multiple 
cavity modes. 



mirror 
x = 




dielectric mirror 
x=q x=l 



FIG. 1: The one-dimensional cavity with moving dielectric membrane at x = ^(t). 



II. THE CLASSICAL LAGRANGIAN AND EQUATIONS OF MOTION 



e{x,q) 



(2) 



We begin by considering a one- dimensional optical cavity of length I formed by two 
perfectly reflecting end mirrors. A movable membrane of rigid uniform dielectric is placed 
inside the cavity. The cavity field is specified by its vector potential A = A{x, t)ez (0 < a; < 
/) under transverse gauge, with the boundary conditions A{x = 0,t) = A{x = l,t) = 0. We 
assume non-birefringent dielectric so that the two polarizations of the field are decoupled, 
and hence it suffices to consider a linearly polarized field. The dielectric is specified by its 
mass m, its center-of-mass coordinate q = q(t)e^ {d/2 < q < I — (i/2), and the dielectric 

constant 

1 + X, q-d/2<x <q + d/2 

1, otherwise 

where d and x are the width and susceptibility of the dielectric, respectively. We have used 
the convention eo = /io = 1 (i-e., c = 1), and assumed non-magnetic dielectric /i = /xq- 
We have also assumed a non- dispersive dielectric so that e does not depend on the field 
frequency. 

The motion of the dielectric affects the electromagnetic fields in the cavity, which in turn 
modifies the radiation pressure on the dielectric. To study the complete dynamics of the 
system, q{t) must be included as a dynamical degree of freedom. The system is specified by 

the Lagrangian 

1 /•' 

L = -mq^ - V{q) + / dxCp, (3) 

2 Jo 



where V{q) is the external mechanical potential on the dielectric, and Cp is the Lagrangian 
(linear) density of the field after eliminating the electronic degrees of freedom of the di- 
electric. To find Cp, we go to an inertial frame S' in which the dielectric membrane is 
instantaneously at rest. Assuming the acceleration of the membrane does not change the 
macroscopic properties of the dielectric, the field Lagrangian density in S' is given by the 
familiar form: Cp = | (eE'^ — B'^) where E' = E'Bz and B' = B'ey are the electric and 
magnetic fields in 5", respectively. Now as the motion of the dielectric in the laboratory 
frame S is perpendicular to both fields q = q{t)ex, we can relate the fields between 5" 
and 5" by the Lorentz transformation: E' = 7 (E + q x B) and B' = 7 (B — q x E), with 
7 = (1 — g^) .In terms of the vector potential, E = —{dtA)ez and B = —{dxA)ey, the 
Lagrangian density Cp in the space between the mirrors reads 

2 / o ^ \ 2 



£, ' ' 



2(1-^2 



(-^^)0-(-^^^)(f)-' 



H^)i^) 



(4) 



This agrees with the earlier work by Barton et al. [28| and by Salamone [29| (in the case 
/i = 1). We see that the Lagrangian density (jlj) appears to be more complicated compared 
with that in the primed frame. This can be understood because in the laboratory frame, a 
dielectric with polarization P moving at velocity q processes a magnetization M = — q x P, 
which must not be neglected in the regime of non-adiabatic dielectric motion. After some 
calculations, it can be shown that the Lagrangian density (j4]) contains the interaction terms 
corresponding to P ■ E + M ■ B. 

In this paper we confine our study for non-relativistic motion q <^ 1, so that Eq. (jl]) is 
approximated by (up to first order in q), 

CF = l[e (dtAf - id^Af] + g (e - 1) (dtA) (d^A) . (5) 

Together with Eq. ([3]), we obtain 

L = Ime - V{q) + I dx{^ [e {dtAf - [d^Af] + g (e - 1) (dtA) (d^A)}, (6) 

which is the non-relativistic Lagrangian of our membrane-field model. 

To justify this Lagrangian, we need to examine whether it consistently generates the 
equations of motion for both the fields and the membrane within the accuracy limited by 
the approximation made in ([6]). First the Euler-Lagrange equation of A{x,t) derived from 



(Q is given by (to first order in q) 

ed^A - dlA + 2g(e - 1) {d:,dtA) + g(e - \)d^A = 0, (7) 

where we have discarded terms involving g^ for consistency, and made use of the relation 
{dt + qdx)^ = 0. From Eq. ([7]), the effects of membrane's motion on the field appear in the e 
terms, the third term with a velocity dependence, and the last term that is proportional to 
the membrane's acceleration. If q = 0, then Eq. ([7]) is simply the wave equation obtained by 
transforming the wave equation ed^,A' — d^iA' = in S" frame to S frame up to first order 
in q 
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32| . The acceleration dependent term therefore acts like a source term in the wave 
equation. However, for an oscillating membrane with a mechanical frequency Vt and field 
frequency w, the ratio of the acceleration dependent term to the velocity dependent term in 
Eq. ([7]) is of the order Vt/u, which is much smaller than one. 

Next, the Euler- Lagrange equation of motion of the membrane based on ([6]) is given by 

q-d/2 

(8) 



dV{q) I X 
mq = ;- h -- 



dq 2 1 + X 



dA^ ^ 



dx 



q+d/2 

The second term on the right-hand side of Eq. ([8]) corresponds to a radiation pressure force 
from the field. Such a force term is consistent with that obtained from the Lorentz force 



density f = (S^P') x B' appearing in S' J33|]. This can be shown by using P' = xE' and the 
wave equation in dielectric rest frame, then a straightforward transform on the force to the 
laboratory frame in the non-relativistic limit would yield the same radiation force expression 
in Eq. ([H]), apart from a term that is about q/c times smaller. Note that our Lagrangian ([H]) 
gives a wave equation that is accurate up to 0{q), but the accuracy is lower by one order of 
q for the membrane's equation of motion. This is because of the partial derivative d/dq in 
the membrane's Euler-Lagrange equation. 

III. THE HAMILTONIAN AND QUANTIZATION 

The Hamiltonian associated with L is defined by 

H (n. A, p,q)=pq+ f dx [n {dtA) - L {A, d^A, q, q)] , (9) 

Jo 



where p and Il{x,t) are canonical momenta conjugate to q and A{x,t), respectively, 

dx{e-l){dtA){d,A) (10) 



P 



n 



dL 

^— = mq + 
dq JO 

""^^ eidtA) + qie-l)id,A). 



(11) 



d{dtA) 

We see that the dielectric canonical momentum p is not equal to its kinetic momentum mq 
for non-zero fields. The explicit expression of the Hamiltonian (Q now reads 

"i rTT2 



H=^{p + Af + V{q) + ^Jdx ^ + {d,A) 



with A given by 



A = - 



dx 



1 



'0 V e 
and m' is identified as a 'renormalized mass' defined by 



U{d,A), 



(12) 



(13) 



m 



m 



m Jo e 



dA\- 
dx J 



(14) 



The form of Hamiltonian ( TT2l) is similar to the minimal coupling Hamiltonian in electro- 
dynamics with A somehow playing the role of the vector potential in the kinetic energy 
term. 

At this point we would like to comment on the renormalized mass m' defined in Eq. 
dH]) ^3^- First, it might look peculiar that the renormalized mass m' depends only on the 
magnetic field energy inside the membrane, but we point out that this is an artifact due 
to the truncation of the Lagrangian (jlj) up to q. If we retain q^ terms in (jlj), it can be 
shown that m' appears to depend on the electric field energy as well. Second, in quantum 
theory, the vacuum field energy would make the integral in Eq. (1141) divergent if all the 
field frequencies are counted. In practice, however, a physical dielectric membrane must 
become transparent (i.e. e — ?■ 1) at high field frequencies, so there is only a finite range of 
field frequencies contributing. It is useful to estimate the order of magnitude if the field 
frequencies are counted up to u^ = 10^^ Hz in the ultra-violet range. We then find that for a 
/ = 1 cm cavity, the vacuum contribution is of the order 10~^^ kg, many orders of magnitude 
lighter than a pico-gram membrane in typical optomechanical setup. If the cavity is filled 
with photon excitations in a single mode, then a similar consideration shows that the photon 
number has to be as high as 10^^ for the mass correction to be comparable to the mass of 




x = l 



FIG. 2: A sketch of a few mode functions. Evidently the mode functions depends on the position 
of the dielectric. 



the dielectric. Hence the mass correction can be safely neglected as long as we restrict our 
dielectric model to optical field frequencies and a sufficiently massive membrane. From now 
on, we will take m' ~ m, and the Hamiltonian reads 



* = i(^+^)' + ''(') + ^ 



dx 



TT2 

— + {d.Af 



(15) 



To quantize the system, we promote the dynamical variables g,p, 74(x), n(x) into opera- 
tors by postulating the commutation relations [g, A{x)] = [q, tl{x)] = \p, A{x)] = [p, tl{x)] = 
0) [QjP] = "^^j [^(3^)) n(x')] = ih6{x — x'). The quantum Hamiltonian takes the same expres- 
sion as flTSl) . but with A defined in flT3|) symmetrized as 



A(g) = - I dx 
'0 



2e 



U d,A]+ d,A]U 



(16) 



A. Instantaneous normal-mode projection 



The field operators can be projected onto any set of complete orthonormal modes. For 
instance, we may use the set of mode functions {ipk{x, go)} defined by 

d'^ipk{x,qo) 



dx"^ 



+ e(x, go)wA:(?o)v5fe(a;, go) = 



(17) 



with a vanishing boundary condition at x = and x = I. Here go is a reference position (c- 
number), say, an equihbrium position of the dielectric. The orthonormahty relation between 
these mode functions is written as: /^ dxe{x, qo)ipk{x, qo)ipj{x, go) = ^kj- Note that we have 
explicitly labelled the mode functions and frequencies as ipk{x,qo) and Uk{qo) to emphasize 
their dependence on go, i.e., ifk{x, go) would be the normal-mode of the field if the membrane 
had been fixed at x = go. 

By substituting A{x) = J2k Qk<-Pk{x, go) and n(x) = Xl^ Pk(^{x, qo)<-Pk{x, go), the Hamilto- 
nian reads 

r n 2 

H 



1 
2m 



k.j J \ k,j k / 



where 



Qk 


= / dxe{x 
Jo 


,qo)(pkix,qo)A{x) 


(19) 


Pk 


= / dx(pkix,qo)Uix) 
Jo 


(20) 


Vkjiq) 


= / dxe'^{x,q)e^{x,qo)ipkix,qo)ipj{x,qo) 
Jo 


(21) 


ikM) 


= — dx 

Jo 


'e(x,g) - r 

. ^{x,q) . 


e(x,go)<^fc(x,go) ^^ 


(22) 



In this way the standard commutation relations: [g, Qk\ = [g, Pk\ = [Pi Qk\ = [Pi Pk\ = 0, 
[Qj.Pk] = ih6kj are preserved. Since v9fc(x,go) is not an instantaneous normal-mode of the 
cavity, the field part of the Hamiltonian (TT8l) is not diagonalized. 

To reveal the physical picture it is always desirable to cast the field part of the Hamiltonian 
in a diagonal basis. This can be achieved by performing a unitary transformation H' = 
T^HTi (where Ti is defined in Appendix R|) : 

2 



H' 



1 
2m 



1 



P+J29kAQA +viq) + -J2\Pk+^liq)Ql 



k,j 



with 



(jm ■ (x Q) 

9kj{q) = - I dx ^ e{x, q) — ^— ^ h [e(x, g) - 1] 



d^j{x,q) 



Vk{Xiq) 



(23) 



(24) 



iQ I (>'y ux 

which satisfies guj = —Qjki so that gkk = 0. The field part of H' is now diagonalized with 



g-dependent frequencies ujk{q)- In particular, the transformed field operators are given by 



tIA{x)Ti = ^ irlQkTij ipk{x,qo) = ^QkVk{x,q) 

k k 

T/n(a;)Ti = ^\TlPkTijt{x,qQ)ipk{.x,qQ) = ^Pkt{x,q)ipk{.x,q) 



(25) 
(26) 

k k 

according to Appendix [XI ipk{x,q) is defined in Eq. flTTl) but with go replaced by q. Noting 
that q is the position operator of the membrane, ipk{x,q) corresponds to an instantaneous 
normal-mode of the cavity. Therefore Qk and Pk become the expansion of transformed A{x) 
and n(x) using the instantaneous normal-modes respectively: 

Qk 



dx 



dx 



tIA{x)Ti e(x,g)v3fc(x,g) 
T/n(x)Ti (fk{x,q)- 



(27) 
(28) 



It should be noted that Qk and Pk are the same operators defined in flT9l) and fl20|) and hence 
they are independent of q. Therefore, the g-dependence of Ti(g), e(x, g) and ipk{x,q) must 
have 'cancelled out' each other in (1271) and (l28l). 



ak{q) = 


1 


UJk{q)Qk + iPk 
uJk{q)Qk - iPk 


V 2nuk{q) 


4(^) = 


1 1 


V 2nujQ) 



B. Generalized Fock spaces 

To represent the quantum state of the system, we introduce the g-dependent annihilation 
and creation operators for each cavity field mode: 

(29) 
(30) 

which satisfy the commutation relation [afc(g), a](g)] = Skj- Since afc(g) depends on g, for 
each position of the dielectric we have a set of Fock states associated with that position. 
These states can be labelled as \{nk},q), where {rik} = {ni,n2,n3, . . .} denotes the occu- 
pation number of each photon mode. \{nk},q) is a simultaneous eigenstate of the photon- 
number operator a|,(g)afc(g) and the position operator g i.e. a|,(g)afc(g)|{^fc}, g) = nk\{nk}, q) 
and q\{nk}, q) = g|{^fc}, g)- Such a set of eigenstates are orthonormal and complete, so that 
any quantum state of the whole system |\E') can be expanded in the basis of these eigenstates 
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I.e. 

/l-d/2 
C(K},g)|{n,},g)rfg, (31) 

where C({nfc},g) is the probabihty amphtude. 

With the help of the g-dependent annihilation and creation operators, the Hamiltonian 
( 123|) becomes 

H' = ^^^ + V{q) + Y^ h^q) Ua, + ^ , (32) 

where we have used a shorthand a^ = afc(g) for convenience, and 

The vacuum field energy appearing in ( l32l) leads to the Casimir force on the dielectric (e.g. 



r(g) = - Y 5Z 9k,j{q) W ^ (afcOj - aja] + a]ak - ala^ j . (33) 



see 



35l | for calculations of a setup similar to our case). We may replace the vacuum energy 
by the Casimir potential energy. However, since the Casimir energy is feebly small compared 
with V{q) in typical optomechanical experiments, its effect should be negligible. 

C. A unitary transformation 

The canonical momentum operator p in the Hamiltonian fl32l) differs from the kinetic 
momentum mq due to T{q). We may apply a unitary transformation to the Hamiltonian 
to make the two momenta coincide. This procedure is analogous to the case of atom-field 
interaction, where one transforms the minimal- coupling Hamiltonian into the electric-dipole 
form under electric-dipole approximation. The transformation is defined with the unitary 
operator 



i / dq'T{q') 

"- Jqo 



(34) 



T2(g) = exp 
The transformed Hamiltonian H = T2 if' T2 reads (Appendix [B]) 

H = h u{q) + ^ hukalak + Y^ [^i^^ (alaj + a]akj + ^i]^ [ala] + flfcOj j , (35) 



k,j 



11 



where 



«(^) = V{q) + ^ E ^ (^'^(^) + 2dt^(^)) (36) 

k 



^if (<?) = -^V^^j^ 



'2\kj 1- 7 ^ \k^lj ± i^kj — ^jk) 



(37) 



Afej(g) — fkj + —^ 2_^ fkhfhj + Try 2^ fkhfhifij + ■ ■ ■ (38) 



2! ^^•'^"•'"^ 3! 

ft h,l 



g 



fkM) = / dq'guM')- (39) 

The Hamiltonian (!35|) is the main result of this paper, which determines the couphng 
strengths ^L- once the mode functions are known. We stress that it is apphcable to general 
motion of the membrane, since no assumption of the motion (except g ^ 1) have been made. 
It should be noted that the Hamiltonian fl55]) contains a photon-number non-conserving part 
proportional to {a\aj + dkO-j), which is responsible for dynamical Casimir effect [36j. 

We remark that the transformation has modified the mode function associated with ak, 
as the transformed cavity field operators (related to that in flTSj) by the combined unitary 
transform defined via T = T1T2) read 

k 

T^n{x)T = -tJ2 V ^ («'^ - 4) e(a;, q)Mx, q). (41) 

A; 

In other words, the mode functions have been changed to (pk{,x,q) = ipk{x,q) + 
J2j ^jk{(j)Vj{^^ q) instead of ipk{x, q). We show in Appendix [Bl that {(pk{,x^ q)} indeed forms 
an orthonormal complete set of mode functions. 

D. Linear approximation 



The Hamiltonian H in (13511 exhibits nonlinear feature in the field-membrane coupling. 
In most practical situations, the potential V{q) bounds the dielectric membrane about an 
equilibrium position go such that x^ = ^ — ^o is small compared with the field wavelengths 
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concerned. Therefore we can perform the expansion 

duk 



OJk{q) ~ ^kO + Xr^ 



akiq) ~ ako + 



Xr. 



dq 

duk 



9=90, 



2a;fco V 9q 



''fcO 



9=90, 



\ ^ - (0) 



(42) 

(43) 
(44) 



where Uko = cokiq = Qo) and gL = gkjiq = qo) are the frequency and couphng constant 



associated with the equihbrium position respectively, and 

, l^koQk + iPk 



OfcO 



(45) 



2huJko 
is the annihilation operator linearized in uJk{q). The linearized 0^0 commutes with both q 

and p. The linearized H reads 

^2 



H = ^ h u{xm) + ^ hukoalQako + x^Fq, 



(46) 



where Fq is the normal-ordered radiation pressure force 



h 



It. ^— \ 

9 2^ 



dujk 



dq 



9=90 . 



r . /WfcO (0) 

Okj + uJkO\ -—Qkj 



'^feoS'O "*" '^koO-jO + Ct/fcO^O + Ctjo'^fcO ) • (47) 



This agrees with the work in Ref. J3l|, and the corrections to the coupling term are of the 
order xi,. 



IV. SOME PHYSICAL CONSEQUENCES OF THE HAMILTONIAN 

Current experimental and theoretical study on optomechanical systems are mainly in the 
regime of single-mode and adiabatic approximations. While our formulation does not require 
these approximations, we show in this section how our Hamiltonian model can be reduced 
to simpler forms when these approximations are applicable. In addition, we indicate some 
possible modifications due to the interaction with multiple cavity modes. 

A. Frequency shift in single-mode limit 

Single-mode approximation is applicable when the cavity field is dominantly contributed 
by a single mode k, and photon excitations in other modes are negligible. As gkk = 0, the 
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single- mode consideration immediately gives F ?^ from (133|) . and so the Hamiltonian (132 
becomes 



H' ^ h u{q) + huk{q)al{q)ak{q). 



(48) 



By linearizing in x„ 



q — qo and employing rotating- wave approximation (RWA), we have 

duk 



-f^' ~ ^ hu{q) + huJkoalQako + x.^ ih 



O'ko'^kO, 



(49) 



9=90, 



where Uko and a^o are the linearized frequency and annihilation operator defined in the 
previous section. Equivalently, (H9|) can also be obtained from (H6|) and (H7I) under single- 
mode approximation together and neglecting the counter-rotating terms. 

The above Hamiltonians fHS]) and (H^ are what one would expect when the dielectric 
motion is adiabatically slow. In obtaining f HS]) we have simply 'dropped out' the j ^ k 
terms in f l32|) . hence scattering processes between photon modes are neglected. However, it 
should be noted that while only the field mode k dominates, a photon in that mode can 
make (virtual) transitions to other modes then back to the k mode. Such a process induces 
a shift in Uk{q), which is the leading order correction to ( l48l) . To determine the correction, 
we work on our full multi-mode Hamiltonian ( 135|) and examine the Heisenberg equation of 
motion for aj{q), 



aj = -tujjaj + 7 I g, , r^^^j 



1 duj + t 1 duj . 
o-i + c^i TT^q 



' ojj dq 



•' ujj dq 



'E feV 



,{+) 



A+) 



«^+(4~^+4^ 



We have used the equation of motion q = p/m and the relation 



daj 
dq 



1 duj t 

-CL ■ 

2uj dq ^ 



(50) 



(51) 



which follows from ( 129|) and ( 130|) . The terms in (1501) that contain a- and aj are fast- 
rotating, and can be neglected in the spirit of rotating- wave approximation (RWA). Provided 
further that the oscillation frequency of the membrane Q is low compared with the frequency 
difference \uk — oojl, we may adiabatically eliminate aj {j ^ k) in the equation of motion of 

[^kj +ijk ) 



ak 



uk + 2^1:^ 



k, 



1 {^k + 2eit^) - {uj, + 2e 



(+) 



ak- 



(52) 
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Hence as a leading order non-adiabatic correction, the frequency Uk in (HHj) should be mod- 
ified into Uk + Afc, where 

2 



(A+) , .(+) 
.(+) v;^ V ■'' ■^''^ 






We emphasize that this frequency shift is caused by the interaction with other field modes, 
as is evident in the expression of (!53l) . Hence the shift is essentially a multi-mode effect, 
even though it is calculated under the single- mode limit. 

If the dielectric membrane is bounded by a harmonic potential u{xm) ~ '"^^^2;^/2, with 
the membrane equilibrium position go at an extremum of ujk{q) (for example, at go = V^), 
the linearized single-mode Hamiltonian under RWA reads 

^2 



n 1 

H ^ h -mQ^xl^ + hukoaloako 



+hxl 



^kO I V^ (0)2 i^kO — ^jo) i^kO + ^jo) 



T+E4 



AuJko^jo 



aloOfco (54) 



where cj^o = i.'^'^^k/dq^) . The last term of ( J54l) comes from the Taylor expansion of 
Afc, which has a leading order of x^. It is intriguing to note that the last term can be 
viewed in two ways: it can be regarded as a shift in field frequency proportional to x^ 



(hence proportional to phonon number under RWA |20|), as well as a shift in the mechanical 
frequency Vt that is proportional to the photon number. 

Eq. (!54|) shows that besides w^^q, the frequency shift is also contributed by the presence of 
other field modes. To compare the effect between these two contributions, we have performed 
numerical calculations based on parameters on the experimental setup of [20!] where / = 0.06 
m, n = y/l + x = 2.2, rf = 50 nm. With go = 1/2, a mode with uouq ~ 1.77 x 10^^ Hz 
(corresponding to A ~ 1064 nm) is present with w^q ~ —3.68 x 10^ Hz nm~^. From our 
calculation, modes with non-zero g)^- contribute about 0.22 Hz nm~^ to the sum in (1M|1 . 
While the contribution from each mode is feebly small, we found that their contributions 
are approximately the same over a wide spectral interval. Since these contributing modes 
have a frequency spacing of about 3 x 10^° Hz, so if we include modes of 10^^ — 10^^ Hz in 
our model, the correction would sum to ~ 0.7 x 10^ Hz nm~^, which is only a few times 
smaller than the contribution from w^q. However, we must point out that our result depends 
on the number of modes included in the model. For this reason, the use of our calculation 
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is limited to an order-of-inagnitude estimate. A more detailed analysis should include the 
dispersion effect of the dielectric, so as to address the spectral interval unambiguously. 

B. Resonant mode transitions 

Resonant mode transitions occur when the mechanical frequency of the membrane Q is 
comparable with the frequency spacing of neighboring field modes. As an illustrative exam- 
ple, let us consider the case where q = qo = //2 and a harmonic potential u{xm) ~ mQ'^x'^/2. 
If the index of refraction of the membrane is sufficiently high, the eigen-frequencies of cav- 
ity modes distribute as doublets. Now suppose the mechanical frequency Q is close to the 
frequency difference of two cavity modes in a doublet, say, ki and k2, then the two modes 
can be resonantly coupled. Neglecting other non-resonant cavity modes, we approximate 
the Hamiltonian (146|) as 

H = hVt [h^h-\ — I +h{ ujia[ai + bj2a}2a2 ) + hr] {h + h^^ i a\a2 + a^ai j , (55) 

where ui {(JJ2) and ai (02) are the linearized frequency and annihilation operator of the 
ki {k2) mode, respectively, b is the annihilation operator of the dielectric motion, 77 is the 
coupling frequency defined by 






Note that due to symmetry of the system, 

duJki (q) 



dq 



d(^kM) 



Q=<10 



dq 



0, (57) 



9=90 



hence the usual radiation pressure term fH^ is absent, and the field-membrane coupling 



to first order in Xm describes the scattering between the two field modes (as in 26|]). In 
particular, the case oi Q = U2 — (jJi (assuming (X'2 > Wi) corresponds to a resonance at which 
the mode coupling can be resonantly enhanced. If a rotating-wave approximation is made, 
the interaction terms in Eq. (155 p would take the same form as that appears in parametric 
down conversion in nonlinear optics. 

As a remark we note that the field frequency spacing can be reduced by increasing the 
cavity length or tuning the refractive index and thickness of the membrane. Recent study 



2l| also shows that avoided crossings of transverse field modes (due to broken symmetry 

16 



of the cavity along its lateral dimensions) can provide a frequency spacing of the order ~ 1 
MHz, which is comparable to the mechanical frequency achievable in current optomechanical 
experiments 37|]. Therefore it would be possible for the membrane frequency to match the 



field frequency spacing. 



V. CONCLUSION 



To conclude, we have presented a non-relativistic Lagrangian and a Hamiltonian for a 
one- dimensional coupled membrane-field system in the g <^ 1 regime. The classical equa- 
tions of motion of both the field and membrane are obtained within the accuracy limited 
by the approximation made in the Lagrangian. Our Hamiltonian ( 135|) should capture op- 
tomechanical processes not described by the single-mode adiabatic approximation and linear 
approximation. For example, we have indicated that the presence of multiple cavity-modes 
can modify the single-mode Hamiltonian [see Eq. ( 15^ ]. and possibly give rise to parametric 
down conversion type interaction if the membrane frequency matches the frequency spacing 
of cavity modes. With the explicit form of interaction strengths between the membrane 
and various cavity modes in fl35|) . one can further study quantum dynamics resulted from 
optomechanical coupling in non-adiabatic regimes. 

There are interesting subtle features that emerge in developing the Hamiltonian model. 
First, the velocity-dependent coupling in the Lagrangian (^, which is necessary to recover 
the leading radiation pressure force term in ([8]), causes the membrane to have a canonical 
momentum different from its kinetic momentum. The unitary transformation T2 can elim- 
inate the difference, but it turns out that T2 also changes the field operators accordingly. 
In this paper we have shown that the transformation T2 on the fields can be interpreted as 
a modification of mode functions. Another subtle feature is the use of instantaneous mode 
functions in diagonalizing the field part of the Hamiltonian fl2^ . The Fock states associated 
with instantaneous modes are defined by the g-dependent photon creation/annihilation op- 
erators. It may not be convenient to perform calculations directly based on such creation 
and annihilation operators, but these operators can always be expanded around the equilib- 
rium position of the membrane in order to obtain the relevant interaction terms. Using our 
approach, such an expansion can be carried out to the first order of Xm (as in Sec. HID), or 
to higher orders as desired. Finally, we should point out that since our model is based on 
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the non-dispersive approximation of the dielectric, a further investigation should incorporate 
the absorptive and dispersive properties. This requires an extension of the current theory 
of quantized field in dispersive media |38| to a moving media. 
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Appendix A: The Ti{q) transformation 



The unitary operator Ti{q) is defined as 
Ti(g) =exp 
with 



E / dq'CkM'' 



Ti 



k,j -"io 



CkM') 



d<^j{x,q') 

dxe{x, q ) — ipk{x, q ) 

dq' 



(Al) 



(A2) 



Here ipk{x,q') are mode functions defined in Eq. ( TT7I) but with go replaced by q'. The 
integrand (kjW) in (IA1|) is a c-number except at the upper limit q. In other words, the 
integral in the exponential of f lAl|) can be viewed as an anti-derivative of CkjW), evaluated 
at the two end-points. It follows that 



tdT^ 



p' = TipT, =p- ihTi^ =p-j2^ (p,Q^ + Q.p,y^ . 



(A3) 



To consider the transformation on Qk and Pk-, we decompose Ti into a (continuous) product 
of infinitesimal transform 



r.= n 



I? =90 



I ^-^ \PkQj + QjPk 



k,j 



Q-qo 



n ll-Uq'K{q') 



(A4) 



where the product is 'g'-ordered', i.e. the leftmost term is associated with q' = qo, the 
rightmost with q' = q. One can check that ( JAip and ( JA4I) are equivalent by expanding the 
product of flA4p into sums of integrals. Each infinitesimal transform acts on Qk and P^ as 



h 



I + '-dq'K{q') ]Qk(l- ^dq'Kiq') ) = Qk + dq'J^ Ckjiq')Qj 



h 



h 



I + ^dq'K{q') ]Pkil- ^dq'Kiq') ) = P,-dq'J2 CAq')Pj 



h 



(A5) 
(A6) 
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with correction only in second order infinitesimals. Hence the infinitesimal transform mod- 
ifies A{x) and 11 (x) by 

I + ^dq'K{q')^ \y2QkM^,q'yj [l - '-dq' K{q')^ 

^ Qk'^'ki.x, q') + dq' ^ QjCkM')Vk{.x, q') 

k k,j 

,d(pkix,q') 



J2Qk 



ipk{x,q') + dq' 



dq' 



y^^Qk'-Pk{x,q' + dq') 



(AT) 



and 



I+'-dq'K{q')^ (y2Pke{x,q')M^,q'yj (^I - ^dq'Kiq')^ 
^ Pke{x, q')(pkix, q') - dq' ^ Pje{x, q')Cjkiq')^kix, q') 

k k,j 



i:p^ 



e{x, q)ipk{x, q) + dq — (e(x, q')ipj{x, q')) 



^ Pk^ix, q' + dq')ipk{x, q' + dq'). 



(A8) 



We have used the completeness relation 



^ e(x, q')(pk{x, q')^k{x\ q') = 6{x, x'), 



(A9) 



and integration by parts 

dx'e{x', q')^^^l^^^,{x', q') = - f dx'^ [e{x', q')M^', q')] ^,{x', q'). (AlO) 
^q Jo ^1 

It follows that combining all the infinitesimal transform readily gives f l25|) and ( !26|) . which, 

together with f lA3|) . lead to ( 123|) . Note that we need not symmetrize PkQj in (!23|) because 

gkj = —Qjki and hence Qkk = 0. This assertion can be proven by performing integration by 

parts on (!24l) . and noting that the dielectric constant e is fixed in its rest frame, so that 

dt'e = 7 (9f + qdx) e = 7g {dq + d^,) e = 0, leading to dqe = -d^e. 
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Qk -- 


= Qk + Y^ \kjQj 

j 


Pk -- 





Appendix B: The T2{q) transformation 

With the unitary transform T2(q) defined in (1M|) . the canonical momentum operator p 
transforms as 

T|pT2 = p-Y, 9kM)PkQj = p - r (Bl) 

SO that p becomes the kinetic momentum mq in the new basis. The transform field variables 
Qk = T^QkT2 and Pk = T^PkT2 reads 

(B2) 
(B3) 

with Xkj defined in f l38|) . Hence the fields in the new frame (i.e. T^A{x)T and T^Il{x)T) 
can be expanded in terms of Qk and Pk by ( HOj) and ( HTj) respectively. 

Next, we show that ip{x,q) is orthonormal and complete. Using the fact [Qj,Pk] = 
[QjjPk] = i^Skj for unitary transformation, together with (lB2p and (lB3p . we have the 
identity 

Afcj + Xjk + 2_^ ^ki^ji = 0. (B4) 

Furthermore, by examining the form of Xkj with fkj = —fjk, it can be shown that 

/ ^ XkiXji = 2_^ MkXij. (B5) 

These two properties of Xjk readily lead to 

dxe{x,q)(pk{x,q)ipj{x,q) = / dxe{x,q)ipk{x,q)(pj{x,q) = 6kj, (B6) 

^ e{x, q)(pk{,x, q)(pk{,x\ g) = ^ e(x, q)ipk{,x, q)(pkix', q) = 6{x, x'), (B7) 

k k 

hence {(pk{x^ q)} indeed forms an orthonormal complete set of mode functions. 
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